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ABSTRACT 

The radial distribution of luminous {Lx > 10^^ erg s"') X-ray point sources in the bulge of 
M31 is investigated using archival Chandra observations. We find a significant increase of the 
specific frequency of X-ray sources, per unit stellar mass, within 1 arcmin from the centre of 
the galaxy. The radial distribution of surplus sources in this region follows the pi law, suggest- 
ing that they are low-mass X-ray binaries formed dynamically in the dense inner bulge. We 
investigate dynamical formation of LMXBs, paying particular attention to the high velocity 
regime characteristic for galactic bulges, which has not been explored previously. Our calcu- 
lations suggest that the majority of the surplus sources are formed in tidal captures of black 
holes by main sequence stars of low mass, M» < 0.3 - OAMq, with some contribution of NS 
systems of same type. Due to the small size of the accretion discs a fraction of such systems 
may be persistent X-ray sources. Some of sources may be ultra-compact X-ray binaries with 
helium star/white dwarf companions. We also predict a large number of faint transients, both 
NS and BH systems, within ~ 1 arcmin from the M31 galactic centre. Finally, we consider 
the population of dynamically formed binaries in Galactic globular clusters, emphasizing the 
differences between these two types of stellar environments. 

Key words: galaxies: individual: M31 - X-rays: binaries - X-rays: galaxies 



1 INTRODUCTION 

It is a well known fact that the ratio of the number of low mass X- 
ray binaries (LMXBs) to stellar mass is ~ two orders of mag nitude 
highe r in globular clusters (GCs) than in the Galactic disc dClarkl 
1 19751) . With the advent of Chandra and XMM-Newton, studies of 
X-ray point sources in external galaxies have become possible, and 
have shown that also there globular clusters are especially abun- 
dant in LMXBs. This is attributed to dynamical processes, through 
which LMXBs are formed in two-body encounters. Due to the pi 
dependence on the stellar density such encounters are frequent in 
globular clus ters and are ne gligible in the field. Currently, th ere are 
13 LM XBs jLiu et alj|2001.) in the 150 globular clusters ( iHarrisI 
1 19961) known in the Galaxy. 

In the central parts of massive galaxies, the stellar densities 
can reach values similar to the densities in less luminous GCs. 
Except for the very inner parts, these densities are still an order 
of magnitude smaller than the densities found in the most lumi- 
nous GCs, where the LMXBs are preferentially found. However, 
the large volume compensates for the smaller density and LMXBs 
can be formed near the galactic centres in two-body encounters in 
non-negligible numbers. Whereas dynamical interactions in globu- 
lar clusters have been intensively investigated, the parameter range 
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typical of galactic centres remains unexplored. Due to an order of 
magnitude higher stellar velocities, the character of the dynamical 
interactions and relative importance of difi^erent formation channels 
in the galactic centres differ from those in globular clusters. 

Due to the large stellar mass contained in the central region 
of a galaxy, a number of "primordial" LMXBs formed through the 
standard evolutionary path exist there too. Although these can not 
be easily distinguished from the binaries formed in two-body en- 
counters, an argument of the specific LMXB frequency (per unit 
stellar mass) can be employed, in the manner similar to the one that 
led to the discovery of dynamical formation of binaries in globular 
clusters. The volume density of the primor dial LMXBs fo llows the 
distribution of the stellar mass in a galaxy t Gilfanovll200 4') whereas 
the spatial distribution of th e dynamically form ed binaries is ex- 
pected to obey the pj/v law jpabian et alj|l97^ . Hence the latter 
should be expected to be much more concentrated towards the cen- 
tre of the host galaxy and reveal themselves as a population of "sur- 
plus" sources near its centre. 

M31 is the closest galaxy with a bulge density large enough to 
host a number of LMXBs formed through dynamical interactions. 
At a distance of 780 kpc jStanek & Gamavichiri998l : iMacri et al.l 
2001) X-ray sources can be easily resolved with Chandra, even 
near the centre of the galaxy. It has been studied extensively with 
Chandra and we use these observations to explore the radial distri- 
bution of bright X-ray point sources in the bulge. The results of this 
study are presented in the Section |2] where it is demonstrated that 
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Figure 1. The radial distribution of the X-ray point sources in M3 1 , exclud- 
ing globular cluster sources and subtracting expected level of CXB source 
density (shown by the dashed hne). The histogram shows the distribution 
of primordial LMXB sub-population as traced by the stellar mass distribu- 
tion. The normalization of the latter is from the best fit to the data outside 1 



the specific frequency of X-rays sources increases sharply inside 
» 1 arcmin. The possible nature of surplus sources is discussed in 
section [3] The details of dynamical formation of binaries in dense 
stellar environments and dependence on the stellar velocity disper- 
sion are considered in the section |4] The results of this section are 
applied to the inner bulge of M3 1 and to the Galactic globular clus- 
ters in sections |5] Our conclusions are presented in the section[6l 



2 RADIAL DISTRIBUTION OF THE X-RAY POINT 
SOURCES 

With the currently available Chandra data it is possible to study the 
spatial distribution of the X-ray point sources in the bulge, without 
being affected by incompleteness, down to the limiting lumionos- 
ity of 10^^^ erg s"'. We restrict our analysis out to a distance of 12 
arcmin from the centre and combine 26 ACIS observations with 
telescope pointings within the central 10 arcmin region of the M31 
bulge for a total exposure time of 201 ks. Details of the data anal- 
ysis, the source lists and the luminosity functions of various su b- 
populations in the bulge are presented in IVoss & Gilfaiiovl j2007h . 

We model the radial distribution of the X-ray sources by 
a superposition o f prim ordial LMXBs and CXB sources, as in 
IVoss&Gilfanovl ( l2006h . The spatial distribution of the former is 
assumed to follow the stell ar mass distribu tion of the galaxy, as 
traced by the K-band light ( lGilfanovll2004h^ We used the K-band 
image of M31 provided by 2MASS LGA Jjarret et al.ll2003h . The 
distribution of CXB sources is assumed to be flat on the angular 
scales of interest. Before proceeding with the fit we removed the 
contribution the from sources other than primordial LMXBs and 
background galaxies. Firstly, we removed 4 identified foreground 
sources, 1 supernova remnant and one extended source. Secondly, 
we excluded X-ray binaries associated with globular clusters, as 
their origin and spatial distribution are different from the "field" 
LMXBs. 

Amo ng our X-ray sourc es 13 are coincident with confirmed 
GCs from lGalleti et all j2004h and 8 with GC candidates. We esti- 
mated the number of random matches by displacing the sources by 
10 arcsec in 4 directions. We found an average of 0.25 coincidences 
with confirmed GCs and 1.0 with GC candidat es. It is well known 
that th e inner parts of M3 1 are depleted of GCs. lBarmbv & Huchral 
( l200lh estimate that 70 per cent of the GCs within 5 arcmin from 



the centre of M3 1 have been detected, leaving ~ 1 6 GCs undetected. 
As only a fraction, ~ 1/5, of the GCs in M31 contain LMXBs the 
expected contribution of LMXBs from undetected GCs is ~3. Due 
to selection effects, the majority of the undetected GCs are of low 
luminosity (absolute visual magnitudes V > -7), and as LMXBs 
are preferentially found in high luminosity GCs the actual number 
of LMXBs in undetected GCs is expected to be < 1. A large frac- 
tion of the GC candidates are not real globular clusters. However, 
an association with an X-ray source raises the probability of the GC 
candidates actually being GCs considerably. We therefore remove 
these sources from our source list too, noting that all conclusions 
of this paper remain unchanged if the analysis is performed with a 
source list in which these sources are included. 

We fit the relative normalizations of the LMXBs and CXBs, 
using the maximum likelihood (ML) test. The best fit is given by a 
model, in which the normalization of CXBs is zero, meaning that 
all the sources are LMXBs. As an alternative, we performed a.x^- 
fit on the binned data, with > 15 sources in each bin, and obtained 
the same result. The probability that the data can be a realization 
of the model is 0.06 using the Kolmogorov-Smirnov (KS) test, and 
6 • lO""* for the;^^ test. The KS test is less sensitive to deviations 
at the end of a distribution, and therefore the result of the;^'^-test 
is more restrictive. We conclude that the LMXB-l-CXB model is 
rejected. 

The visual examination of the data (Fig. [T]l suggests that the 
reason for the rejection of the model is an overdensity of sources in 
the inner 1 arcmin region of M31. Motivated by this we did ax^ fit 
of the same model to the distribution outside 1 arcmin. The best fit 
value of the normalization of the CXB component gives the total 
number of 26±9 sources CXB sources in the entire r < 12 arcmin. 
This value is consistent with the expect ation of 29 backg round 
galaxies, estimated from th e soft band of ,Moretti et al] ( 12003) , us- 
ing the method described in IVoss & Gilfanov! l l2006l) . We therefore 
fix the normalization of the CXB component at the value corre- 
sponding to 29 sources. This gives a total number of the LMXBs of 
64+7 in the entire r < 12 arcmin image. The^^-value is 2.63 for 3 
degrees of freedom. The best fit model is shown in Fig. [T] together 
with the observed distribution. 

Using the best-fit model it is possible to investigate the dis- 
tribution of sources in the inner 1 arcmin and quantify the excess 
in the surface density of the sources. The total number of sources 
detected in the the r < 60 arcsec region is 29. The extrapolation 
of the best fit model into this region predicts 8.4±0.9 sources, and 
therefore the number of surplus sources is 20.6+5.5. The error in 
the latter estimate accounts for the Poissonian uncertainty in the to- 
tal number of sources inside 60 arcsec and for the uncertainty of the 
best fit model normalization. As it is obvious from Fig[T] the con- 
trast between the observed number of sources and that predicted 
from the K-band light distribution increases towards the centre of 
the galaxy. Inside r < 15 arcsec, for example, 9 sources are de- 
tected with only 1.1 sources predicted. The formal probability of 
such an excess to happen due to statistical fluctuation is ~ 3 • 10"*, 
assuming Poissonian distribution. 



3 ORIGIN OF THE SURPLUS BINARIES 

Non-uniform extinction, peaking at the centre of M3 1 , could cause 
the distribution of the A'-band light to deviate from the distribu- 
tion of stellar mass. This possibility can be excluded, however, as 
the extinction towards the centre of M31 is low, Av=0.24 mag and 
A/=0.14 mag lHai]| 19961) , which extrapolated to the A'-band gives 
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Figure 2. The radial distribution of the "surplus" X-ray sources, computed 
as a difference between the data and best fit model shown in figure [T] The 
histograms show distributions of the confirmed globular clusters (thick grey 
line) and globular cluster candidates (thin grey fine). 

Aic=0.03 teinnev & Merrifielj Il998 l). Moreover, a non-uniform 
extinction distribution would also cause non-uniformity in the ap- 
parent colours of the stellar population, which is not observed 
jWalterbos & Kennicutt|[T987l) . 

The surplus sources can be high-mass X-ray binaries associ- 
ated with star formation in the inner bulge of M3 1 . We derive upper 
limits for the star formation rate and the n umber of HMXBs from 
the Hff and FIR luminosities reported by lOevereux et al.l jl994h . 
The combined Ha luminosity of the nuclear region and from dif- 
fuse emission inside the star forming ring (which lies at a radius 
~50 arcmin, i.e. much larger than the region analysed in this paper) 
is 4.3- 10^' erg s ~' (transformed to th e distance of 780 kpc used in 
this paper). From lGrimm et al] ( l2003h we find that this corresponds 
to 3.2 HMXBs with a luminosity above 10^*^ erg s"'. The FIR lu- 
minosity in this region is 5.25-10* Lq, which corresponds to 5.9 
HMXBs with a luminosity above 10^* erg s"'. 

It should be stressed out, that the region these luminosities 
refer to is almost 20 times larger than the region analysed in this 
paper. Moreover it is very likely that the main part of the Ho- and 
FIR emission is not associated with star formation, as the number of 
O-type sta rs is a factor of ~ 200 lower than what would be expected 
otherwise jPevereux et al. To conclude, the HMXB nature 

of the sources in the inner bulge can be excluded. 

The surplus sources could have been created i n globular clus- 
ters that remain undetected. In the catalogue of Gall eti et al. Il l2004h 
there are 64 confirmed GCs hosting 13 LMXBs in the region we 
analysed. The fraction of GCs containing X-ray sources is there- 
fore 0.2. This number is larger than what is found in other galaxies 
jSarazin et al ] |2003l; iMaccarone e7ai]|2003 l). due to the better sen- 
sitivity of our stud y, but consistent with the results for the inner 
parts of galaxies in lKim et al.l ( l2006h . Attributing the x 20 surplus 
sources to undetected GCs would therefore indicate that ~ 100 un- 
observed GCs exist in the inner 1 arcmin region of M3 1 . This is 
much larger than allowed by the completeness level of the present 
studies of GC population in M31, con sistent with only a few un de- 
tected globular clusters in this region tearmbv & HuchrdbOOlh . 

In a related scenario, the surplus sources could have been cre- 
ated in globular clusters at larger distances from the centre of M3 1 . 
Due to the mass segregation the globular clusters drift towards the 



centre of the bulge, where they are destroyed, leaving behind rem- 
nant LMXBs. This scenario has been motivated by Fig[2]where the 
radial distribution of surplus sources is compared with that of the 
globular clusters. Indeed, for a GC of mass 10' Mq the mass seg- 
regation timescale is ~ 10' yr at a radius of 5 arcmin and ~ 10'" 
yr at a radius of 12 arcmin ( ISpitzeJ[l969l) . Assuming that a neu- 
tron star turns accreted matter into radiation with an efficiency of 
~ 0.2, the lifetime of an LMXB is < lO'mrf/Lsv yr, where mj is 
the mass of the donor star at the onset of mass transfer expressed 
in solar masses, and L37 is the average luminosity of an LMXB in 
units of lO'" erg s"'. Taking into account that on average ~ 1/5 of 
GCs in M31 contain LMXBs, a destruction rate of ~ 100 globular 
clusters per Gyr is required to explain » 20 sources observed near 
the centre of M3 1 . This number is comparable to the total number 
of GCs within the entire region analysed in this study, and is ~ 30 
per cent of the total number of GCs in M3 1 . As GCs are not contin- 
ually formed in large numbers in M31, the globular cluster system 
of this galaxy will not be able to sustain such a destruction rate and, 
consequently, the population of X-ray sources observed in the inner 
bulge, for any significant length of time. 

Finally, the stellar density in the central part of the M3 1 bulge, 
~ lO'* Mg/pc', is high enough that LMXBs can be formed through 
dynamical interactions in the same manner as in globular clusters. 
In the following sections we investigate this possibility, and apply 
it to the population of X-ray sources in the inner bulge of M3 1 and 
in globular clusters in the Milky Way. 

4 DYNAMICAL INTERACTIONS IN DENSE STELLAR 
ENVIRONMENTS 

There are three main channels of dynamical LMXB formation op- 
erating in dense stellar environment^} 

(i) In a tidal capture of a neutron star (NS) by a non-degenerate 
single star, a close passage of the two stars induces oscillations in 
the non-degenerate star, and the energy for this is taken from the 
orbital energy. If the energy of the oscillations exceeds the origi- 
nally positive orbit al energy, the stars are captured in a bound orbit 
l lFabianetalll975l) . 

(ii) A collision between an NS and an evolved single star on 
the subgiant or red giant branch (RGB) or the asymptotic giant 
branch (AGB) can lead to the formation of an X-ray binary, in 
which the donor star is a white/brown dwarf or a helium star, de- 
pending on the evolutionary s tage of the evolved star before the 
collision jlvanova et ahlbOOS b*). In the case of a white dwarf donor 
an ultra compact X-ray binary is formed. In this scenario, orbital 
energy is transferred to the envelope of the evolved star, which is 
expelled, leaving the NS and the core of the evolved star in a bound 
orbit (Verbunt 1987). 

(iii) In an exchange reaction, an NS exchanges place with a 
star in a pre -existing binary during a close binary-single encounter 
l lHillslll976l) . 

In the context of LMXB formation in globular clusters in the 
Milky Way, the attention has been initially drawn to the tidal cap- 
tures jFabian et alj|l97^ . while the potential importance of the 
two other m echanisms has been realized few years later faillsl 
ll976l;IVerbuntiil987.) . The estimates of the LMXB production rates. 



We begin with consideration of encounters involving a neutron star. These 
results are extended to formation of black hole binaries in the section |43] 
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which followed, revealed that each of the channels could give sig- 
nificant contributio n to the population o f LMXBs found in th e 
Galactic GCs (e.g. IVerbunt & h"u3 1 19871: lOavies & Be"n3 Il995h . 
Since then large amount of work has been done to understand 
physics of stellar encounters in detail, explore their parameter space 
and derive accurate mathematical prescrip tions for the crossec- 
tions and rates ( e.g.[Press & Teukols kvl 1977 : Lee & Ostriker 1986; 
McMillan et al. 1990; Rasio & Shapiro 1991; Davi es et alj|1992l ; 
Sigu rdsson & P hinney 1993; Heggi e & Hut 19 93). However, there 
have been surprisingly few studies making specific predictions of 
numbers of dynamically formed binaries which could be directly 
compared with their observed population in Galactic GCs. With the 
exception of a few studies considering a handful of individual GCs 
( lDavies&Ben3ll995l ; llvanova et alT2005h) . the rates are usualy 
computed for a set of the representative values of parameters (such 
as stellar density, velocity dispersion etc.) and then extrapolated to 
the entire Galactic globular cluster system (e.g. Verbunt 1987). An 
often used assu mption is also that the number of LMXB s is propor- 
tional to pl/cr„ jVerbunt & Hulll987l ; |Poolev et alJlOOsh . Although 
acceptable as an initial approximation, it is too crude to perform a 
quantitative comparison of the theory with observations. Another 
major limitation of the most of these investigations is that the sub- 
sequent evolution of the newly formed binary, into and through the 
X-ray active phase, is ignored. Due to lack of the attention and ef- 
fort in this direction, it is currently unclear if any of the channels 
strongly dominates over the others in real globular clusters. Even 
less understood is the operation of these processes in the environ- 
ment of the galactic centers. 

It is the goal of this paper to fill these gaps. In particular, a spe- 
cial attention will be paid to the following aspects of the problem, 
which have often been ignored in the previous publications on this 
subject: 

(i) calculation of the encounter crossections and rates in the high 
velocity regime typical for the galactic centers and investigations of 
their velocity dependence 

(ii) critical review and comparison of the stellar environments 
(present day mass function, metallicity, abundance of compact ob- 
jects etc) in globular clusters and galactic centers and investigations 
of their impact on the overall LMXB production rates via different 
LMXB formation channels. Reasonably accurate calculation of the 
encounter rates in the Galactic globular clusters and inner bulge of 
M31, based on their structural properties. 

(iii) account for evolution of the newly formed binary before 
and during the X-ray active phase and estimates of the expected 
numbers of LMXBs based on the derived encounter rates 



4.1 General considerations 

In the following we compute cross-sections and rates of the three 
formation channels, consider their dependence on the velocity dis- 
persion of the stars and discuss various factors affecting their effi- 
ciency in the high velocity regime. 

In section[5]we use these results to calculate theoretical forma- 
tion rates and numbers of observable LMXBs in the bulge of M31 
and in the Galactic GCs. 

Each of the processes depends on the rate of encounters be- 
tween two types of objects, which in a unit volume is given by 
niiioy, where 



where rti and ni are the number densities and Mi and M2 are the 
masses of object type 1 and 2, respectively, (t(v,(,;) is the cross- 
section of the encounter, and F(v„/) is the distribution of relative 
velocities at infinity. Assuming that the velocity distributions of the 
two kinds of objects are both Maxwellian and have the same three- 
dimensional velocity dispersion cr,,, the distribution of relative ve- 
locities is given by 



F(Vrd)dVrd ■■ 
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Due to the effect of gravitational focusing, the cross-section for two 
objects to pass within a distance D of each other is given by 

Cr(Vrel) =nD^\l + —J (3) 

' rt'l ' 

In most cases, the gravitational focusing (the second term in the 
brackets) dominates. Only for very fast encounters (large D and/or 
v,(,;) is Dv^,^i > 2G(Mi + M2). If D is independent on the relative 
velocity, 7 oc lv,.d for slow encounters, and y oc p^v„/ for the fast 
ones. 

Several remarks are in place, concerning the subsequent evo- 
lution of the newly created binary system with a compact object. 

Capture of a neutron star in a bound orbit with a companion 
will lead to formation of an X-ray binary provided that the compan- 
ion star will fill its Roche lobe and mass transfer will commense 
within a reasonable time, shorter than ~ 5 - 10 Gyr. If the ini- 
tial binary separation is too large for this to occur immediately af- 
ter the capture, it can be decreased in the course of evolution of 
the binary. There are 3 main mechanisms, which affect the orbital 
separation: (i) magnetic braking, (ii) gravitational braking and (iii) 
binary-single interactions. The former two are familiar from the 
standard theory of the binary evolution (see Ivan den Heuveill 19921 
for a review). For the companion mass in the 0.3 - I.OMq range 
they will bring the system in to the Roche lobe contact within 5 
Gyr if the initial orbital separation does not exceed ~ 3.0 - 1 .QRq 
and ~ 2.5 - 2>.QRq respectively. The braking mechanism due to in- 
teraction of the binary with single "field" stars is specific for high 
stellar density environments. Its properties are briefly summarized 
below. 

When considering evolution of a binary due to binary-single 
interactions it is conventional to divide the binaries into soft and 
hard, depending on the ratio of their binding energy to the kinetic 
energy of the single star at infinity (Heggie 1975). Soft binaries 
have relatively wide orbits, and interactions with single stars tend 
to widen the orbit further or to ionize the binary. Hard binaries, 
on th e contrary, are on ave rage hardened by encounters with single 
stars l lHut&Bahcallll983t) . The effect of this is that over time most 
binaries with a separation above a critical value are disrupted, while 
the compact ones become more compact. The boundary between 
soft and hard binaries depends on the stellar velocity dispersion 
and the mass ratios and ranges from a ~ 300 - 1000i?Q in a typical 
globular cluster to a ~ few Rq in the high velocity environment of 
the M31 bulge. Due to a linear dependence of the crossection on 
the binary separation the coUisional braking is mostly important at 
wide binaries, where magnetic braking and gravitational radiation, 
decreasing as inverse power of the binary separation, are inefficient. 

The initial orbital separation in the binaries produced through 
tidal captures and collisions with RGB/AGB stars is small and a 
large fraction of this systems will start mass transfer (i.e. become 
X-ray sources) soon after their formation. Only a small fraction of 
them (<20 per cent in GCs and < 2 per cent in M31) will ex- 
perience close encounters with single stars significantly affecting 
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thery semimajor axis, therefore binary-single interactions are not 
an important factor in their evolution. Binaries created through ex- 
changes, on the contrary, typically have wider orbits, and the effects 
of encounters can be important. 

If the initial binary separation is large and the braking mech- 
anisms are insufficient to start Roche-lobe overflow, this can oc- 
cur when the donor star evolves off the main sequence, as a re- 
sult of its expansion during the giant phase. In these systems the 
accretion disc is large and X-ray emission from vicinity of the 
NS is insufficient for the irradiation t o keep the entire disc ion- 
ized, and they are therefore transient tang et al.|[T997h . Further- 
more, mass transfer can only occur while the donor is on the RGB, 
which makes the lifetime of such systems short. While they may 
account for bright sou rces detected in massive elliptical galaxies 
jPiro & BildstenI 20021) . they are too rare to make a significant con- 
tribution to our sample. The NSs in these systems are spun up to 
become millisecond pulsars, and in the Gal actic G Cs a large num- 
ber of these have been observed iLorime 3 l2005h . After the outer 
layers of the giant star have been ejected, a binary consisting of a 
white dwarf and an NS remains. However for the vast majority of 
the systems the binary separation is too large for mass transfer to 
begin, and they will therefore not become observable in X-rays. 



(MS) or horizontal branch (HB) stars and with evolved stars on 
the RGB or AGB. Due to the different structure of the stars, the 
outcome of a collision with a NS is different. Simulations indicate 
that coll isions between an NS and an MS star tend to destroy the 
MS star jPavies et alJl99A . We expect the same to happen to stars 
on th e HB, as their structure is similar to that of MS stars (Dormanl 
I1992I) . Collisions of this type are not interesting from the point of 
view of formation of X-ray binaries. As the envelope of stars on the 
RGB or AGB is less strongly bound to the core, a collision with an 
NS can lead to the envelope being expelled. The outcome is a short 
period binary consisting of the core of the evolved star and the NS. 
If the evolved star had a degenerate core, an ultra-compa ct X-ray 
binary (UCXB) with a white dwarf donor w ill be formed ( VerbunJ 
ll987l : llvanova et al.l2005bl : lLombardi et al.l2006i) . In either case an 
X-ray binary can be created. 

The maximum value of distance at periastron Rcoii, for which 
significant amounts of material can be exchanged in an encounter 
of an NS with a non-degenerate star (with radius R,) is between 
R, and the orbit al separation at which the star fills its Roche-lobe 
l lEggletonlll983ll : 



(6) 



4.2 Single-single encounters 

The formation rates of LMXBs due to tidal captures and stellar 
collisions can be found by integrating the encounter rate, given by 
equation [T] over the relevant parts of the mass f unction of star s 
f(M). We assume that the latter follows the form o f lKrouDjj200lh . 
a broken powerlaw with slope 1.3 from O.I-O.SMq and slope 2.3 
above 0.5 Mq, and is normalized according to 



Mr: 



f(M)dM = 1.0 



(4) 



where M„,„-, is the maximum initial mass of stars that have not yet 
evolved to become stellar remnants at present, and M^o is the lower 
cut-off mass. The number density of stars is then given by n, = 
^j^, where p, is the stellar mass density. We assume the mass of 
all neutron stars to be IAMq, and that they are formed from stars 
with initial mass in the range 8.0 - 30.0Mo. The number density 
of these can then be expressed as = fns-^^ = fiu^,, where 

/„., = f^'^'^dM (for M,„ = O.IMe and M„,„, = l.OMo fns = 

0.0068). We define the rate integrated with the mass function 



yf(M)dM 



(5) 



where integration is performed in the relevant initial mass range 
(see below) and y is from equation [T] With this definition n,n„,r 
gives the rate of encounters in s"' cm"'. For the calculation of F it 
is necessary to know the current radius of a star R(M) as a function 
of its initial mass, as well as its evolutionary stage, which is used to 
define the mass limits of the integral. These i nformations we take 
from stellar isochrones of lGirardi et alj ( l2002h . 



4.2.1 Collisions 

We define an encounter between two stars as a collision if the stars 
come so close that considerable amounts of material is exchanged 
between them, and hydrodynamical effects become important. For 
the collisions with NSs, that are relevant for dynamical formation 
of LMXBs, we distinguish between collisions with main sequence 



where q = M,/Mi^s- For encounters with NSs, this separation 
ranges from ~ 5.4/?, to ~ 2.8/?, for stars with masses in the 
0.1 - I.OMq range. SPH Simulations of stellar enc ounters have 
show n that for M, ^ IMq, the value of R^^u/R, ~ 1.8 jPavies et al.l 
Il992h . which we adopt as the standard value. The value of RconlR, 
given, the encounter rate can be calculated from equations 1 113 1 and 
[5] (but see below regarding the choice of integration limits in equa- 
tion[T). 

When considering collisions between NSs and evolved stars, 
it is important to note that the envelope of a star on the RGB/AGB 
is loosely bound to the core, and the orbital energy of the two stars 
at infinity can be comparable to the binding energy of the enve- 
lope. It is therefore possible that the envelope is expelled without 
carrying off^ enough energy to leave a bound system, or that in the 
high-velocity encounters, the duration of the interaction is too short 
for enough energy to be transferred from the NS to the envelope. 
While simulations indicate that the RG envelope is promptly dis- 
rupt ed, instead of ejected through a common enve lope (CE) evolu- 
tion ( iRasio & Shapiroll99lLlLombardi et al.l2006l), the energy con- 



siderations are similar Adopting the formalism o f lWebbinkI ( 119841) 
and de Kool ( 1990), we assume that the envelope of the RG is 
ejected, and that energy for this (the binding energy of the enve- 
lope Etind) is taken from the orbital energy of the two stars which is 
therefore changed by AEorb- Allowing energy to be lost, e.g. as ra- 
diation or as some of the envelope is ejected with a velocity higher 
than the escape speed, an efficiency parameter a^e is defined, so 
that Ebiijii = a„A/iorfc. The binding energy of the RG is given by 



1.0 GM„„.M 



where M and R is the mass and radius of the 



RG, respectively, and Mp„,, is the mass of the envelope of the RG. 
/i is a factor that relates the simplified equation to a precise inte- 
gral of the gravitational bin ding energy and intern al energy in the 
envelope of the RG, see e.g. lPewi & TaurisI ( l200(]t) . The change in 
orbital energy needed to reach an orbit with a separation Of is given 
by AEorh = - . " vt, + 



2 M+M„s "rel ' 2af 

the RG, M„s is the mass of the NS and v„; is the relative velocity of 
the two stars at infinity. 

For a given encounter velocity, we can now find the final sep- 
aration of the binary by solving 
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<J, [km/s] 

Figure 3. The rate y of NS-RGB/AGB encounters that lead to the formation 
of a binary as a function of the 3D velocity dispersion for 4 values of o-^/l 
(equation|7). The values are 0.01,0.1,0.5,1.0 from the top to the bottom. The 
shaded area corresponds to the range of velocity dispersions in the central 
parts of GCs, wheras the vertical dashed line corresponds to the velocity 
dispersion in the bulge of M3 1 . 



GM,„,M 



1 MM„ 



GM„,,M„ 



2M + M„. 
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(7) 



When we calculate the rate of collisions with evolved stars, ycon, the 
integral over velocities (equation[T]( is only carried out for velocities 
Vrei < V„,ax.c, defined such that the final separation a/ < SRo (this 
choice of the maximum separation ensures that the gravitational 
braking will be efficient on the formed binary) 



Jcoll 



r 

Jo 



F(v„i, a,,)crcollV,-eldVreI 



(8) 



where a^gu is the coUisional cross-section defined by equation |3] 
with D = Rcoii- In figure [3] we compare the formation rates of 
UCXBs due to collisions between RGB+AGB stars and NSs for 
different values of or^/l. The rates were calculated by integrating 
over all evolut ionary stages of sta rs on the RGB/AGB, using the 
isochrones of jGirardi et al.ll2003) . The details of the calculations 
are given below in section |4!4l It is obvious that the choice of a„A 
is very important in the bulge of M31, with an order of magnitude 
difference between the rates of the highest and the lowest value, 
while the difference is relatively small in GCs. 

SPH simulations indicate that the effectivity in stan- 
dard CE evolution is Qf„ < 0.5 (see iTaam & RickeJ l2006l 
and references therein), and in population synthesis stud- 
ies, valu es of Aa^e in the rang e 0.1-1.0 are most often 
assumed jPortegies Zwart & Yunpels on 199"^; iFrver & WooslevI 
19991 : iHurlev et alj l2002l : IVoss & TaurisI 12001 : iBelczvnski et al.l 
20051) , and this seems to give a good fit to the observed properties 
of the post-CE binary population. 

However there are differences between the standard CE evolu- 
tion and the collisions considered here. In the former, the stars are 
already in a bound orbit, and the energy can therefore be transferred 
to the envelope over a longer period of time during a large number 
of orbital revolutions. On the other hand, in a collision, enough 
energy has to be transferred from the NS to the envelope during 
the first periastron passage, so that the two stars remain bound. 
Especially for high velocity encounters, the CE formalism might 



not be directly applicable, as the timescale for the first passage 
can be so short that the NS passes through the envelope without 
transferring much energy to this. While the low- velocity regime has 
been well investigate d using SPH simulations ( Davies et al. 1993: 
iLombardi et al.ll2006[) . showing that a value of Rcoii/R, = 1.8 is ad- 
equate, no investigation of the high-velocity regime has been per- 
formed. We assume that ace A is in the range 0.1-1.0 (1.0 can be 
considered a very conservative estimate, giving the minimum rate 
of LMXB formation through this process, whereas 0. 1 may be a 
rather optimistic value), with 0.5 being our chosen standard value 
for calculations below. 



4.2.2 Tidal captures 

At periastron distan ces above Rcoii, and up to a few times R,, tidal 
capture can happen. IPress & Teukolskvl (l9T% provided a way of 
calculating the energy absorbed by t he stars assuming n = 3/2 
polytropes and lLee & Ostrikerl ( Il986h extended these calculations 
to other polytropic indices. In this formulation, the energy of os- 
cillations induced in the non-degenerate star of mass M, during an 
encounter with a neutron star is 



AE, 



GM. 



M„, \' y IR^ 



(9) 



where Rp is the distance of closest approach. Only the spherical 
harmonic indices 1 = 2 (quadrupole) and 1 = 3 (octupole) are in- 
cluded, as higher indices give negligible contributions to the energy 
(Lee & Ostriker 1986). The parameter r] is defined as 

jl/2 ^3/2 



1 



M, 



(10) 



M, + M„J \R, 

With the tabulated overlap integrals of iLee & Ostriked (Il986), for 
polytropic in dices n = 3/2 and n = 3, we use the numerical method 
described by IPress & Teukolskvl ( Il977h to calculate Ti(r]). From 
equation ^ it is then possible to calculate the maximum value of 
Rp (we call this Rnj) for which capture will occur, when the mass 
of the star and the relative velocity at infinity is known. We use this 
method for the case where the non-degenerate star is on the MS. A 
polytropic index of n = 3/2 is assumed for the mainly convective 
stars of mass, M, < 0.4 Mq, and n = 3 for stars M, > 0.4 Mq hav- 
ing radiat ive coresH For stars on the red giant branch, we use the 
results o f lMcMillanetaLl ( ll990l) . who calculated R,jj for captures 
of a neutron star by a 0.8 Mg star at various evolutionary stages 
along the RGB. This mass is close to the MS turn-off masses for 
the Galactic globular clusters (GCs), whereas the turn-off mass in 
M31 is higher, ~1 Mq. We use their results directly in both cases. 
As tidal captures by evolved stars give a negligible contribution to 
the overall binary formation rates we did not attempt to perform a 
more accurate computation for the case of M3 1 . We neglect tidal 
captures during the subsequent evolutionary stages. The structure 
of stars on the AGB is similar to those on the RGB, but the time 
spent there is much shorter. The tidal capture rate must therefore be 
lower. The time spent on the HB is also very short, compared to the 
MS lifetime, and although the tidal capture rate may be comparable 
to the RGB rate, it is much smaller than the MS capture rates. 



^ Note that this is different from t he mass limit for a fully convective star, ~ 
0.3Mo fe.g. lSpruit&RitteJI 19831) . used later as the mass at which magnetic 
braking ceases being effective. As properties of deep interiors of the star 
are not importa nt for the tidal c apture process, a higher value is used in the 
latter case (e.g. lRav et al.ll987l) 
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The tidal capture rate 7,;^^; is computed as a rate of encounter 
with the periastron distance R^gu < Rp < /f,,^, i.e. contribution of 
very close encounters resulting in collisions is subtracted: 



Jlidal - 



F(Vi-el, Cr ,,)[cr „j„; - Cr collh'reldVrel 



(11) 



The upper integration limit v„,i,x,t is defined as the velocity at which 
Rcoii = R-tid i-e. the term in square brackets is required to be positive 
inside the integration limits. The cr,Mn; is calculated from equation 
[3] with D = R,id. 

Important for the following evolution of tidal capture binaries 
is the timescale on which the tidally induced oscillations are dissi- 
pated. If this timescale is short (as argued bv lKumar & GoodmanI 
so that a large fraction of the energy is thermalized within 
one orbital revolution of the binary, new oscillations will be in- 
duced at each periastron passage. The binary quickly becomes cir- 
cularized with the final orbital separation (from conservation of an- 
gular momentum): 



(^2G{M.+ M„,)Rp + v,,iRpf 
G(M. + M„,,) 



(12) 



For slow encounters (low v„; or small Rp), a ^ 2Rp. If the dissi- 
pation time scale is too short, the quick conversion of the energy 
(of the order of few per cent of the binding energy of the star) may 
cause the star to expand and lead to a merger. The outcome of the 
dissipation process depends on the thermalization timescale and the 
region of the star where the en ergy is deposited, both factors being 
unkown ( Podsiadlowskil [l99^ . Alternatively, if the dissipation is 
inefficient, coupling of the orbital motion with oscillation can cause 
large fluctuations in the orbital energy, substantially extending the 
circularization process and potentially scattering a fraction of the 
binaries to wider orbits (Kochanek 1992; Mardling 1995). The de- 
tails and the final outcome of this processes are poorly understood. 
In the following we assume that all binaries become circularized 
with the final separation given by equation[T2l 

Equation [TT] defines the total tidal capture rate, irrespective 
of the subsequent evolution of the tidally formed binary. In or- 
der to calculate the rate of encounters, leading to formation of an 
LMXB, one needs to account for the finite braking time scales. 
For this, R,ij, used to calculate anjai in equation[TT] is replaced by 
min(i?„rf, Ri„-ake), where Ri,rake depends on the mass of the star and 
is chosen so that the braking time scale for the tidal capture binary 
is shorter than 5 Gyr. Due to small values of the tidal capture radius 
R,jd, this does not affect the final rates significantly. 



4.3 Binary-single interactions 

The rate of exchange reactions between a binary (Mi , M2, a) and a 
star M3 is: 



j f{M,)dM, j , 



p(q)dq 



, dn 



X I ■y{Mi,MT,M^,a,(r,)—da 
da 



(13) 



where f{M) is the distribution of mass of one of the stars in 
the binary, q is the binary mass ratio and p{q) its probablity 
distribution, dnida is the binary semimajor axis distribution and 
7(Mi, M2, M3,a, cr„/) is the exchange rate of a star M3 into a bi- 
nary (Mi,M2,a) computed from the equation [T] In the context of 
LMXB formation the third star M3 is an NS or a black hole. 

A significant amount of effort has been invested in the past 
decades in studying binary-single interactions and in calculating 



the encounter crossections and rates. The three-body problem in- 
volved in encounters between a binary and a single star can not 
be solved analytically and the computational demand of the nu- 
merical solution has been prohibitive for the studies of large en- 
sembles of binaries based on direct integration. An approach sug- 
gested and successfully implemented in 80-ies was to build large 
libraries of interactions covering interesting range of initial param- 
eters and, based on these libraries, to derive various semi-analytical 
formulae describing t he interaction c r ossection and outcom e (e.g. 
iHut & Bahcail 1 1 9831 : iMikkolal 1 1984 iHeggie & Hu3 Il993l) . This 
provided the basis for computation of the elementary encounter 
rates 7(Mi, M2, M3, a, (T„;). Another ingredient required to com- 
pute the final encounter rates in an ensemble of binary and single 
stars using eq[T3]is the semi-major axes distribution dn/da. This 
distribution has complex time evolution, defined by the counter- 
play of binary (de-)excitation and ionization processes which are 
difficult to take into account analytically, even with the elementary 
crossections and rates given. Furthermore, unlike in collisions with 
red giants and in tidal captures, the binary separation of a typical 
exchange binary is large and one would have to take into consider- 
ation the subsequent evolution of the binary parameters, before the 
Roche lobe contact is achieved and an X-ray sources appears. 

Not surprisingly, the Monte-Carlo methods has been proven to 
be more efficient. In these, the evolution of each individual binary 
is followed through a number of encounters with s ingle stars and, 
in some implementations, with other binaries (e.g. 'Hut & BahcalJ 
1983 ; Heggie&Hut 1993; Sigurdsson & Phinney 1993; Davie] 
I1993I) . Due to computational limitations, these early simulations 
often relied on libraries of interactions and semi-analytical crossec- 
tions, rather than direct integration of the three-body problem for 
each interaction. While making possible to evolve sizable popula- 
tions of binaries, this approach has its deficiencies, as a number of 
distributions functions (semi-major axes, eccentricities etc.) had to 
be replaced by average values or treated in a simplified way (e.g. 
iHut & B^ call 1983). Rapid advance in computing power and nu- 
merical methods in recent years have allowed full simulations, with 
each binary- single star interaction bei ng calculated exactly (e.g. 
iFregeau et alj2004l ; llvanova et al.l2006h . This lifts the assumptions 
and approximations mentioned above, and it is this method that we 
employ in the current st udy. This approach has been implemented 
as early as in 90-ies (e.g. lPortegies Zwart et alJl997h . but the mod- 
erate numbers of binaries achievable then have been insufficient 
to study the formation of LMXBs (see table [TJ. The next level of 
complication (and of computational demand) is a complete time de- 
pendent simulation of an N-body system composed of binary and 
single stars with realistic mass functions. Presently this is becom- 
ing f easible, but is still limited to systems containin g < 100.000 
stars jPortegies Zwart et alj2007l ; lHurlev et alj2007h . 



4.3.1 Monte-Carlo simulations 

In our simulations we follow the evolution of binaries in an environ- 
ment of single stars, with special emphasis on interactions between 
the single stars and the binaries. The outcome gives a reasonable 
indication of the importance of this process, compared to the two 
paths of dynamical formation of LMXBs from single-single en- 
counters discussed above. The simulations are based on the FEW- 
BODY code of Fregeau et al. (2004). FEWBODY numerically in- 
tegrates the orbits of the stars during the interaction, and automati- 
cally classifies and terminates calculations as soon as the outcome 
is unambiguous, which makes it well suited for carrying out large 
sets of binary interactions. 
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This simulation was performed for a limited range of orbital separations and y„ have been corrected for this. 



Table 1. Parameters and results for the three simulations of exchange reactions. The parameters are: the velocity dispersion cr,, in km s" , the number density 
of single stars n, in pc"', the NS fraction in the population of single stars /„, and the number of binaries simulated Ntin - The results are: the total number of 
NS binnaries formed thi'ough exchanges of a neutron star in to a binary (Total) and the numbers of LMXBs formed in the simulations - with MS (MS-NS) and 
RG (RG-NS) donors. The formation rate 7„ calculated from equation [15] (cgs units), is given for LMXBs with MS donors. The other rates can be computed 
by scaling with the numbers of formed systems. The two groups of columns present results with and without account for stellar collisions. 



All binaries and single stars are assumed to be formed at the 
same time, and the simulation of the binaries begins 0.5 Gyr after 
the star formation episode. The mas ses of single st ars are assumed 
to follow the initial mass function of iKroupjHoOll) . Only main se- 
quence stars and neutron stars are included, and the main sequence 
turn-off mass evol ves with the age of the population (as given by 
iGirardi et alj200j) . The number density of single stars is kept con- 
stant during the simulation. 

The initial binaries are drawn randomly from a population 
with properties typical of binary population sy nthesis studies (e.g. 
iDewev & Cordesll 19871 ; IPoIs & Marinuslll994l) . The mass distribu- 
tion of the primary stars {Mp) was chosen to be the same as the 
mass function of single stars, while the mass of the secondary cho- 
sen from a flat mass ratio distribution. The distribution of orbital 
separations a was assumed flat in log a between a minimum sepa- 
ration corresponding to one of the stars filling its Roche-lobe, and 
a maxim um separation of \Q*Rc, con sistent with the distribution 
found bv lDuquennov & Mavoil ( Il99 1[) in this range. The initial ec- 
centricity of the binaries was set to 0. Each binary is evolved for 
15 Gyr in the single star environment, taking into account stel- 
lar evolution and evolution of the binary orbit due to magnetic 
braking and gravitational radiation as well as encounters with sin- 
gle stars (but omitting the more complicated l ate phases of bi- 
nary e volution of more advanced models, e.g . Dewev & Cordes' 
19871: IPols &MarinusI Il994 I Porte gies Zwart & Verbunt 1996; 
Hurley et al.ll2002l ; lBelczvnski et al. 200^^ 



The binaries are evolved in timesteps of maximally 0.01 times 
the average time between encounters with single stars, where an 
encounter is assumed to happen if a star comes closer than 6 
orbital separations. For each timestep binary parameters are ad- 
justed according to gravitational wave emis sion jLandau & Lifshit3 
ll962l : |Peterj|l964h and magnetic braking jRappaport et al.lll983h . 
We assumed the disrupted magnetic braking model, where mag- 
netic braking is ineffective when the MS sta r is totally convective 
jRappaport et alJl983l : [Spruit & Ritterll"983h . This is the case when 
the mass of the star is below ~ 0.3 Mg. The probability of an en- 
counter between the binary and a single star within a timestep of 
length A/ is given by a weighted average over the distribution of 
relative velocities and over the mass function of single stars. 



^n, \ yf(M)dM 
Jm 



(14) 



Here n, is the number density of single stars, y is given by equa- 
tion[T] with the cross-section found from equation [3] with D = 6a, 
where a is the orbital separation of the binary. Random numbers 
are drawn to see whether an encounter occurs. If this is the case, the 
parameters of the encounter are drawn from their respective prob- 



ability distributions. The mass of the single star is drawn from the 
mass function (allowing for neutron stars also). The probability of 
an encounter distance D is proportional to ^ (where cr is given by 
equation[3]l out to the maximum distance of 6 orbital separations of 
the binary. The distribution of encounter velocities Vrei is given by 
equation|2] The encounter is then solved for using the binsingle 
program of FEWBODY. Binary phase and encounter angle is cho- 
sen randomly by FEWBODY from a flat and an isotropic distribu- 
tion, respectively. The simulation of a binary is terminated when 
one of the following occurs: (1) the binary is disrupted, (2) one of 
the stars evolves off the main sequence or (3) Roche-lobe contact 
is reached. If one of the binary components is an NS, possibility 
3 leads to the formation of an LMXB. Possibility 2 also leads to 
Roche-lobe overflow, but as discussed above, such RGB-NS sys- 
tems are shortlived and transient X-ray sources. 

The simulations are performed with several simplifying as- 
sumptions. We discuss the most important of them below. En- 
counters between binaries are ignored. As most wide binaries are 
quickly destroyed, the binary fraction decreases fast and binary- 
binary encounters should only matter at early times. Moreover, 
in most binary-bina ry encounters, two or more of the stars merge 
jFregeau et alj2004}) . Secondly we have negl ected the effec t of tidal 
interaction in the evolution of the binaries ( IZahnlll989al lbl). This 
will tend to lock the rotation of the stars to the orbit and to cir- 
cularize the orbit, thus decrease somewhat the time it takes for a 
system to achieve Roche-lobe contact. The significance of this ef- 
fect is diflicult to e stimate, for its impleme ntation in population 
synthesis codes, see iBelczvnski et al.l ( l2005h . Evolved stars were 
not included in the simulations. For the single star population this 
should not be a problem, as an encounter between an evolved star 
and a binary will probably lead to a merger of some sort due to the 
large radius of the evolved star. The net effect of such encounters 
will most likely be a decreased binary fraction. As for the evolved 
stars in the binaries, they will lead to Roche-lobe overflow. It is 
unlikely that a neutron star can be exchanged into such a system 
without the occurrence of a physical collision. We verified with test 
simulations that close encounters between tight binaries and single 
stars in which one of the stars is evolved in almost all cases lead to 
m erger of two or all thr ee stars, in accordance with the conclusions 
of iFregeau et alj ( l2004h . 



4.3.2 Results of simulations 

We performed three simulations with different velocity dispersions 
and densities, to cover the environment in both M31 and in GCs. 
Parameters and results of the simulations are summarized in table 
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Time [Gyr] 

Figure 4. The evolution of the dynamical formation rate of LMXBs due 
to exchange reactions in a low velocity environment of a globular cluster 
(cr,, = 15 km s"'). The plot shows the number of systems in which Roche- 
lobe contact was reached, per time bin, as a function of time from the start of 
the simulation. Results of simulations without and with account for mergers 
are shown with thick grey and thin black crosses respectively. 



[T] Presented in the table are the numbers of neutron star binaries 
created in the simulations - the total number and the numbers of 
Roche-lobe filling systems. The latter is divided into the following 
two categories: the binaries in which Roche-lobe overflow occurs 
due to evolution of the binary orbit, while the companion star is 
on the MS (MS-NS) and the systems in which the mass transfer 
is initiated due to the evolution of the companion star off the main 
sequence (RG-NS). 

We convert the numbers into 7„ rates, which can be directly 
compared to the single-single interaction rates computed in the pre- 
vious section, using 



Ye, 



N 



(15) 



where and tij^s are the number densities of binaries and neu- 
tron stars, Nbin is the total number of binaries in the simulation and 
Tsim is the simulations time span. Note, that with the above defini- 
tion, Nbin is the total number of binaries simulated, i.e. ni,m has the 
meaning of the primordial volume density of binaries. The rates for 
MS-NS systems are given in the Table [T] The other rates can be 
computed from these by scaling according to numbers of binaries. 
The Monte-Carlo uncertainties can be estimated assuming a Pois- 
sonian distribution for the numbers of binaries. In the M31 simu- 
lations (the largest value of cr,.) we simulated binaries in a limited 
range of separations, as wide binaries are quickly ionized and do 
not contribute to the LMXB production rates. The final value of the 
exchange rates given in Table[T]has been corrected correspondingly. 

In the initial simulations the radii of stars were set to zero, i.e. 
the possibility of stellar collisions was not accounted for. We per- 
formed a second set of simulations, in which two stars with radii 
R\ and R2 were assumed to collide if the distance of closest ap- 
proach was D < 1.8(i?i + Ri)- We assumed that collisions lead to a 
merger and removed from simulations all binaries that experience 
such events. With these assumptions, the rates of LMXB forma- 
tion decrease dramatically, by about an order of magnitude, i.e. in 
~ 90% of binary-single interactions which could potentially lead to 
an exchange of the neutron star into the binary, two or more of the 



stars collide. This result is in agreement with lFregeau et al.] ^2004) 
and demonstrates that mergers are a determining factor in the for- 
mation of exchange LMXBs. 

Figure|4]illustrates the time dependence of the formation rate 
of exchange LMXBs in a low velocity environment of a typical 
globular cluster (simulation with cr,, = 15 km s"'). Shown in the 
figure is the number of binaries per time bin, in which the Roche- 
lobe overflow was initiated during the given time bin, irrespective 
of the time when the neutron star was exchanged into the binary. 
There is an obvious increase with time, due to the fact that in a low 
velocity dispersion environment most MS-NS binaries are created 
with relatively large orbital separations and need to be hardened by 
further collisions in order to become LMXBs. This is in contrast 
to M31 (large velocity dispersion), where the rate is constant, due 
to the fact that almost all exchange LMXBs there are formed from 
binaries with small orbital separations, a < lOi?©- For such bina- 
ries, binary-single interactions are not an important factor in their 
further evolution towards Roche-lobe overflow. Also shown in the 
plot by thin crosses is the result of simulations with account for 
mergers. 

For interactions, where the final binary is harder than the ini- 
tial binary, the binding energy lost is converted to kinetic energy of 
the binary and the single star. The velocity of the binary obtained 
due to thi s effect is often referred to as the dynami cal recoil ve- 
locitv(e.g. lSigurdsson & Phinnevlll993l : lDavieslll993l) . We find that 
in globular clusters, the binaries that end up as LMXBs typically 
undergo several encounters, in which recoil velocities in the range 
~ 30 - 50 km s"' are obtained. The effect of this is that the bi- 
nary is ejected from the core (and sometimes also from the GC) 
to the less dense regions of the GC, where dynamical interactions 
are rare. After a significant time, the binaries will re-enter the core 
due to mass segregation. For a recent discussion of this binary cy- 
cling in and out of the core of a GC using N-body simulations, see 
lHurlevetal] ( l2007l) . As this cycling can lead to significant prolon- 
gation of the binary lifetime before the formation of a LMXB, it 
decreases the formation rate due to this channel. In the simulations 
with physical collisions this effect is smaller, due to the fact that 
many of the encounters that lead to high recoil velocities are also 
the encounters that lead to collisions. We note that the encounter 
cross-sections of LMXBs are so small that only a small fraction of 
them experience significant encounters, and they are therefore re- 
tained in the GC cores. In M31 most LMXBs are formed through 
only one encounter, an exchange reaction in which the orbital sepa- 
ration is decreased significantly, for an already tight binary. The re- 
coil velocities are therefore extremely high, typically 100-600 km 
s"' . Even for the deep potential wells of galactic bulges, such veloc- 
ities can be enough to eject the binaries. However, as will be shown 
in section |5] this channel is not a significant source of LMXBs in 
M31 anyway. 

4.4 Comparison of the rates 

The results of this section are summarized in figure [5] where we 
compare the rates for the three main LMXB formation processes, 
involving neutron stars, as a function of the stellar velocity dis- 
persion. In computing the rates for the tidal capture and collisions 
with evolved stars we assumed an environment (IMF, age metal- 
licity etc.) similar to the bulge of M31, as described in section[5] 
As will be discussed in section [5] the stellar environment in GCs 
is significantly different from that in M31, in particular with re- 
gard to the present day mass function and metallicity, and this is 
taken into account in our final estimates. The goal of this section is 
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to highlight the influence of the velocity dispersion. The exchange 
LMXB rates are from the simulations of the previous subsection, 
without and with account for mergers. It should be noted, that low 
and high velocity parts of these simulations were tailored for GC 
and M31 environment respectively, therefore were performed for 
different values of the stellar density, main sequence cut-off mass 
Mco, age and metallicity. For this reason, although they do correctly 
illustrate the general trend of the exchange rates with the stellar ve- 
locity dispersion, they should not be used to study the exact depen- 
dence. One should also keep in mind, that in order to convert the 
formation rates into the numbers of X-ray sources, the LMXB life- 
time considerations should be taken into account, as discussed in 
the section l4~6l 

Figure [5] illustrates significant velocity dependence of the 
relative importance of different LMXB formation channels and 
suggests that the relative contributions of different subclasses of 
LMXBs should be different in GCs and in the galactic centres. In 
the low velocity environment of a GCs all three processes make 
comparable contributions to the LMXB production rates (but not to 
the numbers of X-ray sources observed in any given time, see be- 
low), with some prevalence of tidal captures by the main sequence 
stars, depending on the exact value of the velocity dispersion. In 
the high velocity environment of a galactic bulge, on the contrary, 
the tidal capture by main sequence stars with M, > 0.3Mq and ex- 
change reactions are unimportant and the LMXB formation rates 
are dominated by the collisions with evolved stars and tidal cap- 
tures by very low mass stars. However, the comparison between 
globular clusters and M3 1 is more complex than comparison of the 
velocity dependent rates, as these environments also differ in other 
properties of the stellar populations, such as the present day mass 
function, metallicity, binary fraction etc. This is considered in detail 
in section|5] 

Finally, the total rate of encounters in volume V can be ob- 
tained as: 



R 



where / = f„s, fu,, fun- Note, that the former two coefficients refer 
to the present day values, while the latter is the primordial binary 
fraction, as clarified in the previous subsection. 



4.5 Black hole encounters 

The discussion above was limited to the formation of LMXBs in 
which the compact object is an NS. Of course, the same processes 
are relevant for black holes (BHs), and these are considered be- 
low. We assume that stars with initial masses in the 30 - IOOMq 
range become BHs, with a canonical mass of 10 Mq. The rates 
of tidal captures and and collisions with evolved stars can then be 
found from the equations of the previous sections, replac ing M„j 
with M j,;, and f„s with fi,i,. With the initial mass function of lKroupal 
( l200llfl fbi, = 0.17 • Note that although there are a 6 times 
fewer BHs than NSs, this is countered by the the gravitational fo- 
cusing term (equation [3]l which makes the encounter cross-section 



^ To use the mass function of lKroup3 jlOOll) to estimate /;,;, //„s is appro- 
priate for the bulge of M31, where mass se gregation is negligible, except 
for the inner few parsec jpreitag et ai]l2006l) . In GCs this is not the case as 
three factors affect the mass function, namely the mass segregation caus- 
ing more massive objects to sink to the core, supernova kicks depleting the 
GCs of NSs and the ejection of B Hs due to encounters with other BHs 
jPortegies Zwart & McMillarj|2002l) . This is discussed in detail in section|5] 
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Figure 5. Comparison of NS LMXB dynamical formation rates as a func- 
tion of the stellar velocity dispersion. Shown are the rates ycr,. for tidal 
captures by main sequence stars > 0.3Mq (thick solid line) and < O.SMq 
(thin solid line), collisions with evolved stars (assuming At] = 0.5, dashed 
line) and exchange reactions (circles with eiTor bars). For the exchange re- 
actions, results are shown for simulations without (open symbols) and with 
(filled symbols) account for collisions. Error bars are only shown if they are 
larger than the circles. The hatched area shows the velocity dispersion range 
typical for globular clusters, the dotted line for the inner bulge of M3 1 . 
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Figure 6. Comparison of the LMXB formation rates in encounters with 
neuti'on stai's (thin lines) and black holes (thick lines). Shown are tidal cap- 
tures by main sequence stars of mass > 0.3Mo (solid lines) and collisons 
with RGB/AGB stars (dashed lines). The normalization of the BH rates have 
been multiplied with fbhl fns- 



~5 times larger for a BH. As for the neutron stars, we assume that 
RcoiiIRt = 1-8. We note, however, that from equation |6] it might 
be expected that may be larger for the very small mass ratios 
considered here. 

When considering various aspects of single-single encounters 
on the basis of the total energy arguments it should be noted that 
the kinetic energy at infinity in the centre-of-masses frame is de- 
fined by the reduced mass ji = MiM2/(Mi + Mi) which, for the 
low mass stars, only depends weakly on the compact object mass. 
Therefore, even if the energy equipartition is inefficient and the NS 
and BH velocity dispersions are comparable, as is the case in the 
M31 bulge, the kinetic energy at infinity, yuv^^,/2, will not be much 
higher in the case of a black hole. For this reason the energy consid- 
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erations in collisions with evolved stars (section 14.2.1 1 eq|7j, will 
not change significantly and the fraction of collisions expected to 
lead to a bound systems for a given velocity dispersion depends 
only weakly on the mass of the compact object. 

For tidal captures, the energy of the oscillations induced in the 
non-degenerate star is roughly oc M^^^^ (see equation |9]l, and the 
capture distance is therefore larger for BHs than for NSs. At low 
velocities, this can enhance the rate by a moderate factor of <2. 
At large velocities, on the contrary, i?„y ~ Rcoii, and even a small 
increase in i?,,^ can drastically increase cr,,^ - crcoii and thereby the 
overall tidal capture rate. As the total energy to be dissipated in the 
interaction depends weakly on the compact object mass, the impact 
of the tidal oscillation on the thermal state of the normal star does 
not become more severe for the black hole. Therefore, if a tidal 
capture is possible at all, it is possible for black holes as well as for 
neutron stars. 

There is one important difference between BHs and NSs, 
namely that the rete ntion factor for black holes in globular clus- 
ters is close to zero jPortegies Zwart & McMillai]l2002i) . For this 
reason, black holes do not contribute to dynamical LMXB forma- 
tion in globular clusters. In M31, however, the black hole fraction 
should be close to the IMF-based estimate given above, due to 
much longer energy equipartition time scale than in globular clus- 
ters. 

In figure [6] we compare the LMXB formation rates in single- 
single encounters involving black holes and neutron stars. Obvi- 
ously, black holes can make sizable contribution to the LMXB for- 
mation rates, especially in the high velocity regime. 



4.6 Numbers of X-ray sources 

In order to convert the encounter rates to the numbers of X-ray 
sources observed at any given moment of time, one needs to con- 
sider the evolution of a binary through the X-ray phase. A definitive 
answer can be obtained from proper population synthesis calcula- 
tions, which is beyond the scope of this paper. In a simpler approach 
one may consider characteristic life times Tx of binaries at different 
phases of its evolution. The number of X-ray sources Nx can be 
then related to the corresponding encounter rate: Nx ~ y^x- Tak- 
ing into account dependence of the Tx on the mass and evolutionary 
status of the donor star and their mass distribution, we obtain an ex- 
pression, similar to the equation[T6]for overall encounter rate: 

A'x « (tm^)' ^ / ^""^^^ "^"^-^^^^ "^^"^^ ^^^^ 

where / is defined as in ea.lf76t and Tx(M) = AM^/M, AM^ = 
Mi — Mf, Mi is the initial mass of the donor star and Mf its final 
mass in the given evolution stage (e.g. for a star with initial mass 
> 0.3Mo, Mf = 0.3Mq - the mass corresponding to the period gap). 
In case of an LMXB powered through the Roche-lobe overflow, the 
M is defined by the orbital braking mechanism and the mass-radius 
relation for the donor star The stability of the mass transfer in the 
accretion disc should be also taken into account in order to identify 
persistent/transient nature of the binary. The integral in ea.lll7t is 
taken over the range of the masses relevant to the given type of 
X-ray binaries. 

Below we examine evolution and characteristic values of M 
of X-ray binaries formed formed via different dynamical processes 
considered in in this paper. We accept th e stan dard prescriptions 
for the m agnetic braking ( Rappaport et al.lll983f) and gravitational 
radiation jLandau & Lifehitzil 9621 : IPetersI 19641) and the transiency 



00 




- 0.2 0.4 0.6 0.8 



Figure 7. Dependence of the mass accretion rate M in a Roche-lobe filling 
system on the mass of the donor star. The curves for a neutron star and a 
black hole binary are shown by thick and thin lines. The calcultaions based 
on 10 and 1 Gyr isochrones give identical result down to ~ O.IMq, be- 
low which the 10 Gyr isochrones give more pronounced spike in M. The 
transiency limits are plotted by respective dashed lines. The method of cal- 
culations and assumptions are described in section l4!6l 

criterium in the form published by iDubus et for irradi- 

ated discs. One should keep in mind that these simple prescip- 
tions predict time averaged quantities but may fail to explain the 
momentary values of luminosity, which may vary significantly 
on the timescales of days-months-years. The dependences of the 
mass accretion rate on the mass of the donor star for NS and 
BH binaries are shown in figure |7] These dependences were com- 
puted b ased on standard formulae for a Roche lobe filling sec- 
ondary jvan den Heuvel 1992) assuming that t he secondary i s in th e 
thermal equilibrium. As was demonstrated bv lStehle et al.l ( Il996h . 
this assumption gives sufficiently accurate results for the main se- 
quence dono r. For the mass-rad ius r elation we used the 10 Gyr 
isochrones o f lBarafreetal1 ( l2003f) and Bara ffe et ai](ll998') for stars 
M < O.IMq and M > O.IM©, respectively. Also shown in figure 
|7] are transiency limits for different types of compact object. The 
NS and BH masses were assumed 1.4 and 10 Mq respectively. The 
spike in M at » 0.07 Mq is caused by the steepening of the mass- 
radius relation just above the hydrogen burning minimum mass, 
as given by the isochrones. This is due to the fact that correlation 
effects between particles becomes important, and the deviations 
from an ideal gas decreases the pressure (I. Baraffe, private com- 
munication). Below the hydrogen burning minimum mass, degen- 
eracy effects dominate and the mass radius relationship becomes 
R ~ M"'''. The spike is less pronounced in the 1 Gyr isochrones 
(shown in figure |7]by thin solid lines) which might be more appro- 
priate for the thermal state of a mass-losing brown dwarf. 

Collisions with evolved stars. In a collision with a red giant, an 
ultra-compact X-ray binary (UCXB) with a He white dwarf donor 
is formed. The white dwarf mass equals app roximately the mass 
of the red giant's core jLombardi et alj|2006h . i.e. is in the 0.1-0.4 
Mq mass range, depending on the evolutionary stage of the red 
giant. The evolution of such a system includes a very fast initial 
stage of very short, ~ minutes, orbital period and very high, super- 
Eddington M. During this period the white dwarf donor is quickly 
reduced to a ~ O.O6M0 after which a more "normal" UCXB with 
forb ~ 10 min and Lx < 10^^ erg s"' emerges, similar to the 
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ones observed in our Galaxy. Overall, such a system will spend 
~ 0.1 Gyr with the luminosity 10^* - 10'* erg s before the white 
dwarf is depleted below » 0.02Mo. Somewhere around this mass 
the s ources will cross the stabi l ity threshold and will bec ome tran- 
sient jPelove & BildsterJl2003l : lBildsten & Delovell2004l) . 

The cores of less evolved, sub-giant stars are not fully degen- 
erate and/or hydrogen-depleted. In this case a collision will result in 
a binary with He or brown dwarf donor, depending on the core mass 
and chemical composition. Such a binary is also driven by gravita- 
tional radiation, but due to the larger radius the period minimum is 
higher, ~ 20 - 30 minutes, and super-Eddington mass transfer does 
therefore not occur. For such systems, a life time of ~ 200 - 300 
Myrs may be expected (N.Ivanova, private communication). In or- 
der to make a crude estimate of their fraction we assume that the 
core of an RGB star becomes fully degenerate, when the central 
density exceeds Pr > (5 - I0)pcrii, where p„,-, is the critical den- 
sity above which electron gas is degenerate (pc™ ~ 2.4- lO^^/j^r^'^ 
g cm We estimated from the Padova stellar tracks that this oc- 
cures at stellar radii of ~ (3 - 5) x RqM/Mq. As discussed in sec- 
tion |4.2. Il given the high stellar velocities in M31, only RGBs with 
rather small radii can elfectively capture a compact object through 
collisions, and we expect that in a large fraction, ~ 50-80 per cent, 
of X-ray sources created through this mechanism the donor star is 
not fully degenerate. In the low velocity environment of globular 
clusters this fraction is smaller, ~ 25 - 40 per cent. 

Tidal captures by main sequence stars with M > 0.3Mq lead 
to formations of "usual" LMXBs, similar to the ones constituting 
the majority of systems with main sequence donors observed in the 
Galaxy. These sources are driven by the magnetic braking and lu- 
minosities of ~ 1036 5-3 8.0 grg/sec and lifetime s of ~ 0.1 - 0.5 Gyr 
should be expected (e.g. Ivan den Heu veil 1993) . Note, that these es- 
timates depend critically on the magnetic braking prescription, the 
weak magneting breaking predic ting up to several times sma ller lu- 
minosities and longer life times jlvanova & Kalogera|[2003) . From 
figure |7] it can be seen that all black hole systems are expected to 
be transient, in agreement with BH binaries statistics in the Milky 
Way. 

Tidal captures by main sequence stars of very low mass, M < 
O.SMq. For these fully convective ver y low mass stars the m agnetic 
braking is believed to be inefficient ( Spruit & Ritteij |l983f) . there- 
fore the accretion is driven by the gravitational radiation. From the 
requirement that the donor star fills its Roche-lobe we have that the 
orbital periods of these systems are in the ~hours range, and that 
gravitational radiation can provide luminosities of ~ lo'*""'*-^ and 
~ 1 036-5-37.0 grg^s fq]- jsjs BH systems res pectively (see also 
e.g. iPodsiadlowski et alj[2003 : lYungelson et alj,200 6). i^From fig- 
ure |7] it can be seen that the systems with M > O.ISM© will be 
transient, these constraints being more severe for the NS binaries. 
Integration of the mass transfer rate gives that the life times during 
the persistent phase are ~ 300 Myr, and that the life times during 
the transient phase are ~ 1 and ~ 4 Gyrs for BH and NS systems 
respectively. 

It is interesting to consider the final stage of evolution of these 
systems, after the donor star is reduced to < 0. IMg, below the 
nuclear burning limit. As these are descendants of very low mass 
stars, whose nuclear time scale is much longer than the cosmolog- 
ical time, they consist mainly of hydrogen and they will become 
brown dwarfs. Given the mass-radius relation for brown dwarfs, 
the mass transfer rate drops quickly when the mass reaches 0.05 
Mq (Fig|7l(, and these systems become transients, similar to some 
of the accreting msec pulsar systems, observed in our Galaxy. 



We note that in the binary systems with very low mass ra- 
tios, q< 0.02, t he circularization radius exceeds the tidal trunca- 
tion radius (e.g. |Paczvnski|[l977l) . It is therefore not entirely clear 
whether the stable mas s transfer is possible, see e.g. discussion in 
lYungelson et al.l ( |2006|) (section 3.3). Such low mass ratios can be 
reached for the most low mass black hole systems. 



5 M31 AND THE MILKY WAY GLOBULAR CLUSTERS 

Below we compute rates of dynamical formation of LMXBs and 
their expected numbers in M31 and in Galactic GCs. For this, we 
need to specify velocity dispersion, initial and present day mass 
functions, age, metallicity and stellar isochrones. These parameters 
are different in GCs and galactic centres. The difference in stellar 
velocities is an important one, as discussed in section l4!4l but sev- 
eral other properties of stellar populations play equally significant 
roles in shaping the population of dynamically formed binaries. 
The factor of prime importance is highly efficient mass segregation 
in GCs. Its two most significant consequences are: 

(i) The present day mass function. Due to efficient mass seg- 
regation, the inner regions of the GCs, where most of the en- 
counters happen, are depleted of low mass st ars, to the degree 
that the mass function is essentially flat (e.g iKing et"aL I I2OO5I : 
Ide Marchi & Parescd 1 19971 : lAlbrow et al] I2OO2I) . This is not the 
case for a galactic bulge, where the mass distribution of main se- 
quence stars is suffi ciently well represented by the Kroupa IMF 
jZoccali et al.ll2000l) . As a result, the tidal captures by very low 
mass stars, dominating the binary formation processes in M31 
(FigjSj, are significantly less important in GCs. 

(ii) Abundance of BHs. GCs are b elieved to be depleted of black 
holes jPortegies Zwart &McMillaij 12002) . due to mass segrega- 
tion and BH-BH encounters (although the observ ation of an ultra- 
luminous X-ray source in a GC in NGC 4472 bv lMaccarone et all 
( 2007), indicates that some BHs may be present in GCs). Therefore 
tidal captures of BHs do not play any roles in globular clusters as 
opposite to the case of M3 1 . Note that in the latter case the role of 
black holes is further enhanced by the velocity dependence of the 
tidal capture cross-section, as discussed in section 1431 and shown 
in figurejS] 

Among other factors, leading to further quantitative differ- 
ences, the following should b e mentioned: (i) Due to supernova 
kicks ( lLvne&Lorirnerlll994l) . large fraction of neutron stars es- 
cape the parent cl uster, with the NS retention fac tor being in the 
~ 0.1 - 0.2 range torukiei|[T99^ : IPfahl et alj[200A . On the other 
hand, the mass segregation of the remaining NSs may increase their 
density near the globular cluster centres, thus compensating for the 
low retention fraction, (ii) Binary fractions are different in globular 
clusters and galactic centres, due to different rates of binary-single 
processes, caused by difference in velocities and stellar densities. 
This is important for exchange rates and is taken into account in 
our simulations automatically, (iii) Finally, different ages and, espe- 
cially metallicities result in different mass-to-light ratios, the main 
sequence turn-off mass and duration of the red giant phase, as dis- 
cussed in section [5!2| 

For these reasons the comparison between globular clusters 
and M3 1 is more complex than comparison of the velocity depen- 
dent rates, shown in figure|5] It is the subject of this section. 
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Figure 8. The radial distribution of "surplus" X-ray sources computed as 
a difference between the data and best fit model in figure [T] The solid line 
shows the projected pi distribution, computed from the original mass model 
of the M3 1 bulge from Rift'eser et al. ( 2005). The dashed line was computed 
from the mass distribution with the circumnuclear stellar disc excluded. 
Both model distributions are normalized to the observed number of surplus 
sources outside 4 arcsec. 



5.1 M31 

For the stellar models we used an ischrone with a metallicity ~ 
1.5 times solar dSalaiedini & Jablonka|[200^ and an age of 12 Gyr 
( lRiclill2004lH This gives a main sequence turn-off mass of 0.9532 
Mq and a mass at the tip of the AGB of 1 .008 1 Mq (initial masses)|f| 
Stars more massive than this have all turned into stellar remnants. 

The velocit y dispersion (3D) was assumed to be constant, 
o-,.=260 km s ' iMcElrovll 19831 : Iwidrow et alj|2003l). The density 



struc ture of M3 1 was constructed using the model of iRiffeser et al 



(2005), based on the Gunn-r band photometry presented by iKent 



(1989). In this model the total iJ-band luminosity of the bulge out 
to a distance of 12 arcmin from the centre of M3 1 is 1 . 1 8 • lO'^Lgj;. 
We normalized the density by requiring the integrated R-band lu- 
minosity over the mass function (giving a mass-to-light ratio of 
(M/L)r=3.27) to match the R-band luminosity of the model, giving 
the bulge a total mass of 3.9-10'°Mo of stars in the 0.1-1.0081 Mq 
range. The projection of this model agrees with the mass distribu- 
tion inferred from the A'-band light, which was used to model the 
LMXB distribution in Voss & Gilfanov (2007). The observed K- 
band luminosity of the region is 4.4-10"'LQji:, and integrating over 
the isochrone, we find a mass-to-light ratio of (M/Qx^O-TG, giving 
a total mass of 3.4-10"'Mo, compatible with the /?-band estimate. 
As a consistency check, we estimate the mass, using the mass-to- 
light ratios of Bell & de Jong (2001). With the bulge colour (B- 
V)=0.95 dWalterbos & Kennicuttlll987h . we find a bulge mass of 
3.75 and 3.73 -IO'^Mq from the /?-band and the A'-band, respec- 
tively. 

In figure [8] the pi profile, integrated over the line of sight, is 
compared to the observed distribution of surplus sources, which 



"* Isochrone file isoc_z030.dat from lGirardi et alj i2002l) . 
^ We note that the precision of these masses is given for identification on 
the published stellar tracks. This should not be confused with the precision 
to which these values are actually know, which is much lower. 



was calculated by subtracting the best-fit model of LMXBs and 
CXBs from the observed radial distribution of X-ray sources (sec- 
tion |2l figure [T). It is obvious that the distributions agree well ev- 
erywhere outsi de ~ 4 arcsec. In the i nnermost 4 arcsec of M3 1 the 
mass model of IRiffeser et al.l ( l2005h features a sharp increase in 
density, absent in the distribution of X-ray sources. This increase is 
due to a stellar disc of high density su rrounding the central super- 
massive black hole (' Bender et al. ?2005). In this paper we do not try 
to model the environment in this region and exclude the disc com- 
ponent. The stellar model used for computation of the encounter 
rates is described by the following distribution: 



Phtdge 

where 



,n-0.4(7.1o,"f +0.61) 



0.254z2 + ^0.2542z*+4(.Y2+y2 + 1.11)^2 



(18) 



(19) 



with at„/je, xo, Jo and zo expressed in arcmin. The inclination of the 
bulge coordinate system is assumed to be 77° , and po = 4.34- 10''Mo 
pc"'' (using our mass to light ratio {M/L)r = 3.27). This gives a 
bulge mass (within 12 arcmin from the centre) of 3.87 ■ IO'^Mq and 



pUv = 4.6- 10"M^pc" 



(20) 



It is now straightforward to calculate the rates of tidal captures 
and collisions. Following the equations of section |4] the rates are 
given by 

^«3i=r i^^f dV ■ U f f{M)rdM (21) 

J bulge \ < M > / J Ml,,,, 

where M/o,,. - M;,,^;, is the initial mass range for the type of stars 
for which the rates are calculated. The rates for different types of 
encounters are summarized in Table |2] For clarity the channels ex- 
pected to lead to the formation of LMXBs are written in bold font. 



5.7.7 Numbers of X-ray sources 

We turn now to the numbers of of dynamically formed X-ray 
sources. As it is obvious from Table [2] (column 4), the number of 
"normal" presistent LMXBs with a neutron star accreting from a 
main sequence companion M, > 03Mq, which constitute the ma- 
jority of the primordial LMXBs, is negligibly small (BH capture 
products with M, > 0.3Mo donors are expected to be transients 
and are discussed below). The two main contributions to the pop- 
ulation of dynamically formed sources come from the tidal cap- 
tures of black holes and neutron stars by very low mass MS stars, 
and from collisions of compact objects with RGB stars (columns 3 
and 5). In computing the numbers of sources from equation [TT] we 
take into account that the the evolutionary timescales of all types of 
dynamically formed X-ray sources are much shorter than the life 
time of the bulge. Therefore the systems formed via tidal capture 
by Mt > 0.3Mq stars will pass through the phase of the very low 
mass companion in the end of their life time, adding to the num- 
bers of persistent and transient sources of this type. Similarly, a 
capture product of, for example, a 0.3Mo star will go through the 
transient phase in the beginning of its X-ray active phase and will 
become a persistent source after the donor star is depleted below 
~ 0.10 - O.ISMq. We thus predict ~ 24 and ~ 5 persistent X-ray 
sources with black holes and neutron stars respectively, accreting 
from the very low mass stars. To this number should be added the 
number of ultra-compact X-ray binaries produced via collisions of 
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Object 


Type 


MS(< 0.3Mo) 


MS(> 0.3Mo) 


RGB 


HB 


AGB 


Exchange 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NS 


Tidal capture 


15.1 (15.7) 


0.8 (0.9) 


0.01 






3.5 


NS 


Collisions 


36.6 


46.6 


6.5-13.3 


5.2 


0.01-0.67 


0.3 


BH 


Tidal capture 


65.3 (65.3) 


14.2 (14.2) 


0.09 






8.8 


BH 


Collisions 


37.7 


38.0 


3.1-5.8 


2.0 


0.00-0.06 


1.0 



Table 2. Formation rates of LMXBs in M31, per Gyr. The entries in bold are expected to lead to the formation of LMXBs. The columns are: (1) The type of 
compact object; (2) the capture process; (3) rate of interactions with MS stars of mass < 0.3Mo, only those initiating Roche-lobe overflow within 5 Gyr are 
included (full number is given in the parenthesis); (4) rate of interactions with MS stars of mass > 0.3Mo, same criteria as (3); (5) rate of interactions with 
stars on the RGB, for collisions only those with af < 5i?o are included for r/A = OA - 1.0; (6) rate of interactions with stars on the HB, tidal captures were 
not calculated; (7) rate of interactions with stars on the AGB, tidal captures were not calculated, same criteria as (5); (8) LMXBs created through exchange 
reactions, estimated from simulations without collisions (in tidal capture rows) and with collisions (in collision rows). 



compact objects with red giants, which is ~ 3. The total number 
of predicted persistent sources is compatible with, albeit somewhat 
larger than the observed number of surplus sources, ~ 21. Given 
the number and magnitude of uncertainties involved in the calcu- 
lations and the simplifications made, we consider this as a good 
agreement. 

Based on the range of the donor masses corresponding to un- 
stable mass transfer (figure |7}, we predict ~ 30 BH and ~ 22 
NS transient sources with very low mass donors M, < 0.3Mq, as 
well as ~ 3 BH transient sources with MS donors > 0.3Mo. Fur- 
thermore, exchange reactions might contribute with a number of 
LMXBs with RGB donor stars, that are also transient, but dura- 
tion of their active phase is restricted by the life time of the red 
giant donor. The number of transients observed at any given mo- 
ment in time depends on their duty cycle. Taking Galactic black 
hole transients with the main sequence donor as an example, one 
could expect a duty cycle of ~ 1/50, giving one bright transient in 
~ 15 years. As for the transients with very low mass donors, one 
can use the accreting msec pulsars as an example of NS systems. 
SAXJ 1808. 4-365 8 has outbursts lasting for ~ 2-3 weeks every ~ 2 
years, and the duty cycle is therefore ~0.03. Assuming crudely that 
it is the same for BH and NS systems, we would expect 1.5 tran- 
sient sources at any given time. The outbursts of accreting msec 
pulsars in our Galaxy are characterized by low peak luminosities, 
log(Lx) 36 - 36.5. Therefore many, if not most, of outbursts 
from these sources will be missed i n a Chandra survey of the type 
reported in IVoss & Gilfajiovl ilOOH) which detects mostly brighter 
transients, with the peak luminosit y of log(Lx) > 36.5. This ex- 
plains whv lVoss& Gilfanov! ( |2007|) have not found any excees in 
the number of transient sources close to the galactic center - the 
fraction of transients detected inside 1 arcmin from the center (5 out 
of 28 in 29 Chandra observations with the time span of ~ 5 years) 
agrees with the fraction of stellar mass contained in this region. On 
the other hand Chandra observations of our Galactic Center, hav- 
ing much better sensitivity, have indeed revealed overabundance of 
faint transients (iMuno et alj2005i) . 

5.2 Globular Clusters 

Due to high efficiency of the mass segregation in globular clusters 
the (retained) neutron stars will be much more centrally concen- 
trated than low mass stars. Assuming that stellar density and veloc- 
ity dispersion are approximately constant over the region occupied 
by the neutron stars, one can approximately write: 

I n„,n,T dV - n,Tc I ihis dV = kN„sn,Tc (22) 
Jv Jv 



where A'„ , is the total number of neutron stars in the globular clus- 
ter under consideration, n, is the central density of stars, is the 
central value of F (equation [5j and k < 1 is a constant accounting 
for inaccuracy of this approximation. Assuming that the distribu- 
tion of normal stars follows the analytic King model (King 196^ , 
and that the NSs are in t hermal equilibrium with the stars at tum - 
ofi' (0.80-0.85 Mo) (as in lLugger et alJll995l ; lGrindlav et alj|2002h . 
we estimated that k ~ 0.2 - 0.3. Total thermal equilibrium is gen- 
erally not reached, the value of k is therefore slightly lower. In the 
following we will use a value of k = 0.2 for all globular clusters. 

We use the catalogue of iHarrisI [m^ for the globular clus- 
ters parameters required to estimate the formation rates of LMXBs 
in the Galactic GCs. Of the 150 GCs included in the catalogue, 
the parameters are missing for 10, and we ignore these. The stel- 
lar populations in t he GCs were modelled using the isochrones o f 
iGirardi et al] i2002h . with an age of 1 1 Gyr jSalaris & Weissl2002h . 
The N„s for each GC was computed as follows. Assuming the initial 
mass function of Kroupa (2001), we used the integrated light of the 
isochrones to compute the present day mass-to-light ratio and from 
the total V-band luminosity of the GCs computed the IMF normal- 
ization. Assuming further that all stars with the initial mass in the 
range of 8 - 30 Mq have become neut ron stars and retention factor 
of 10 per cent (Drukier 1996 ; Pfahl et alj2002h we finally compute 
the present day number of the neutron stars in each globular clus- 
ter, A',,,. On the other hand, we assumed that the present day mass 
function in the GC centers is flat. With this mass function we again 
use the integrated l/-band light of the isochro nes to calcula te n, 
from the V-b and luminosity density p y given in lHarrisI ( Il996i) . For 
the 56 GCs in lPrvor & MevlanI ( [19931) we use their central velocity 
dispersions Vq needed to compute the encounter rates. The remain- 
ong GCs were dealt with as follows. ^From the virial theorem we 
expect that vo ~ Kr^ -s/po, where is the core radius of the GCs, 
and po is the central density; we further assumed that pv oc po- We 
performed the least square fit to the known central velocity disper- 
sions in 56 GCs and found K = 0.18 km s"' and 0.17 km s"' for 
the metal-rich and metal-poor GCs respectively (assuming that 
is in pc and pv in Mq y pc"^). These values have been used to find 
Vo for the remaining 84 GCs. 



5.2.1 Metallicity effects 

In order to study the metallicity dependence of the encounter 
rates, we compute the cumulated rates for two metallicities, 20 per 
cent, and 2 per cent o f the solar va l ue (fi les isocz884.dat and 
isoczO0O4.dat from lGirardi et al.l ( |2002|) ) which are representa- 
tive of the red and blue GC populations, respectively. The results 
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Metallicity 


Type 


MS(< 0.3Mo) 


MS(> 0.3Mo) 


RGB 


HB 


AGB 


Exchange 


0.2 solar 


Tidal capture 


8.5 (10.4) 


29.3 (32.5) 


7.0 






203.4 




Collisions 


5.6 


56.1 


24.3-27.7 


4.6 


0.4-1.1 


15.3 


0.02 solar 


Tidal capture 


5.5 (6.6) 


17.3(18.9) 


2.9 






117.3 




Collisions 


3.6 


31.6 


10.1-11.7 


2.0 


0.3-0.6 


8.8 



Table 3. Total encounter rates for 140 Galactic globular clusters from lHarrisI (1996) for which sufficient structural parameters are known, calculated assuming 
metallicity of 0.2 and 0.02 solar. Entries in bold indicate paths expected to lead to the formation of LMXBs. The rates are given in LMXBs/Gyr These 140 
GCs contain the 13 LMXBs observed in the Galactic GC system. The notation is the same as in table|2] 



Population 


LMXBs observed 


Type 


MS(< 0.3Mo) 


MS(> 0.3Mo) 


RGB 


HB 


AGB 


Exchange 


Red GCs 


8 


Tidal capture 


2.7 (3.1) 


7.9 (8.7) 


1.6 






53.3 






Collisions 


1.8 


18.4 


7.9-9.2 


1.5 


0.1-0.4 


4.0 


Blue GCs 


5 


Tidal capture 


3.8 (4.6) 


12.6(13.8) 


2.2 






86.1 






Collisions 


2.4 


21.3 


6.8-7.8 


1.3 


0.2-0.4 


6.5 



Table 4. Total encounter rates calculated separately for red and the blue Galactic globular cluster subsystems, assuming metallicity of 0.2 and 0.02 solar 
respectively. Bold entries indicate paths expected to lead to formation of LMXBs, and the rates are given as LMXBs/Gyr The notation is the same as in table 

in 



are presented in tableland show a ~L5-2.5 increase in the en- 
counter rates for the higher metallicity case. 

The metallicity dependen ce in our calculations is mainly due 
to two factors. (1) As noted bv lBellazzini et al] ( Il995l ) the radii of 
metal-rich stars are larger, and therefore the rates of tidal captures 
and collisions are higher. Furthermore the duration of the RG phase 
is long er for metal-rich stars. As demonstrated by Maccarone et al. 
( 12004) this effect can maximally lead to an enhancement of the 
cross-sections and rates by < 60 per cent, and most likely ~ 30 
per cent. Our results are consistent with this, showing a ~ 20 per 
cent increase in tidal captures by MS stars (> O.SMq) and ~ 50 per 
cent increase in collisions with RGB/ AGB stars. For exchange re- 
actions the effect is negligible. (2) Theoretical isochrones predict 
that the V-band mass-to-light ratio of the metal-rich population is 
higher than that of the metal-p oor population. As the stellar den- 
sities are given in HarrisI jl996t) in the form of (/-band luminosity 
density, the encounter rate is proportional to pi oc {MIL)-. This 
could result in an additional ~ 60 per cent increase in the rates. It 
is however unclear, whether this is the case for real globular clus- 
ters - observations indicate that the the central mass-to-light ra- 
tio might be independent on the metallicity tMcLaughlinI 120001) . 
This could be due to the fact that the red GCs typically are more 
dynamically evolved (but not older) than the blue ones and there- 
fore have a flatter mass function in their cores (' McClureetalJ 19861 : 
IVesp eiini & Heggie 1997; Piotto & Zoccali 199^ Moreover, these 
structural differences may be the true reason for the observed metal- 
licity dependence of the abundance of dynamically created sources 
in globular clusters as also noted bv lBregman et al. U2006h . 

Thus our calculations do indicate a moderate metallicity de- 
pendence of the encounter rates. It is however insufficient to explain 
observations. Indeed, there are ~3 times as many LMXBs in red 
GCs as in blue GCs of the same size in the Galaxy ('Grind lavll99l : 
iBellazzini et"ai. 1995), where 8 out of 13 LMXBs are observed in 
the red GC system containing 46 out of the total number of 140 
GCs with known metallicities (assuming a divi sion at [Fe/H]=-1) . 
Similar trend is observed in in othe r galaxies jKundu et alj|2002l : 
ISarazin et al]|2003l : iKim et al.ll2006l) . 



5.2.2 Predicted rates and numbers of X-ray sources 

To predict the total rates of LMXBs formation in the Galactic GCs, 
we divide the GCs into two subpopulations depending on metallic- 
ity, red (46 GCs) and blue (94 GCs). The cumulative rates for these 
two subpopulations are then calculated as above, assuming all red 
GCs to have 0.2 solar metallicity and all blue ones to have 0.02 so- 
lar metallicity. The results are given in table|4] As it can be expected 
from figure|5] all three processes give comparable contributions. 

For metal-rich clusters, these rates predict ~ 1.2 X-ray bina- 
ries with the companion mass > 0.3Mq due to tidal captures, with 
an additional 0.5-1.0 such binaries from exchange reactions, ~ 1.5 
UCXBs and ~ 3 fainter LMXBs with very low mass companion. 
Corresponding to ~ 6 sources overall, this is in a good agreement 
with the total number of LMXBs observed in metal-rich clusters 
(8). On the other hand, we do overpredict the numbers of X-rays 
sources in the metal-poor GCs by a factor of ~ 1.5 - although 
our calculations do show the expected metallicity dependence, it 
is compensated by the larger number of metal-poor clusters. Note 
that the number of bright sources with Mj > 0.3Mq main sequence 
companion dependes critically on the rate of magnetic braking. The 
ab ove numbers have been computed with the standard prescription 
of lRappaport et al] (Il983l). The weaker variants of magnetic brak- 
ing (e.g. Ivanova & Taan]||2003l) may give upto a factor ~ 5 - 10 
longer lifetimes and, consequently, larger numbers of LMXBs with 
Mrf > 0.3/Wo donors. This can change the overall numbers for glob- 
ular clusters, but is insignificant factor in the M31 bulge calcula- 
tions, due to negligible contribution of these systems there. 

It is interesting to compare the numbers of ultra-compact sys- 
tems. Considering metal rich clusters only, 2 of the 8 LMXBs 
bital period s 
jistall2003) . 

mined periods and could therefore be either UCXBs or traditional 
LMXBs. The final 2 have orbital periods > 5 h. Thus, there may be 
from 2 to 6 short period systems. We predict ~ 1.5 UCXBs formed 
in the collisions with red giants. In addition, the LMXBs with the 
very low mass donor stars, < Q.\5Mq, for which the predicted 
number is ~ 3, will also have short orbital periods and faint op- 



have mea sured orbital period s < 1 h and are therefore most likely 
UCXBs teenacQuistal |2006|). Of the 6 others 4 have undeter- 
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47Tuc 




oj Cen 






Tidal 


RG-NS 


Tidal 


RG-NS 


Our study 


4 


2 


11 


3 


Davies & Benz (1995) 


3 


1 


14 


2 



Table 5. Comparis on of the LMXB produ ction rates for two Galactic GCs 
with the results of iDavies & Benj jl995l) . When computing our numbers 
we adjusted the parameters of the stellar environment according to the as- 
sumption of lOavies & Benj Im^ . as described in the text. The rates are 
given in units of Gyr"' . 



tical counterparts and may contribute to the observed statistics of 
UCXBs, giving a prediction of ~ 4.5 short period systems in total. 



5.3 Comparison with previous studies 

As we have already emphasized above, there is only a handful of 
studies dedicated the quantitative predictions of the formation rates 
and numbers of LMXBs in real globular clusters. 

The results of our globula r cluster calculations agree with 
the estimates of I Verbung ( l2003h of the relative production rates 
of LMXBs in different Galactic GCs. Later on this method has 
been used to succesfully explain the observ ed distribution of X- 
ray sources ov er Galactic i Poolev et al ] |2003[) as well as for extra- 
Galactic GCs i iSivakofretalJl2007l) . 

One of the most d etailed investigations so far has been done by 
lDavies&Ben3 ( ll995 ') who considered LMXB formation in realis- 
tic models of w Cen and 47 Tuc globular clusters. When compared 
blindly, their results appear to differ from our calculations for the 
same two clusters. However, this is due to different assumptions on 
the stellar environments used in t heir study. The main differences 
are, that in lDavies & Benj ^^3): (i) the NS depletion effect due 
to supernova kicks was not taken into account (i.e. 100% NS reten- 
tion factor has been assumed), (ii) it was assumed that stars more 
massive than 6Mq produce NSs as compared with 8Mq boundary 
used in this paper and (iii) the initial mass function for co Cen was 
much flatter than Kroupa IMF used in this paper. In order to com- 
pare with their results, we modified our calculations to be consis- 
tent with these assumptions and we found good agreement between 
the two studies, as demonstrated by the table |5] We note that there 
also is a number of more subtle differences, not taken into account 
here, which may explain the remaning differences. Finally, consis- 
tent with our results, they find that the production rate due to ex- 
change reactions is similar to the rate from the other two channels. 

Ou r calculations a re also consistent, within a factor of 2, with 
those of llvanova et al.l I^OSb), who studied the rate of collisions 
between RGs and NSs in a small sample of Galactic GCs. 

The only investig ation of the formation of LMXBs in galactic 
bulges is the study bv lLee&Nelsonl ( ll988h of tidal captures near 
the Galactic centre. Extrapolating from calculations of captures by 
stars of 0.5 Mq only, they found that while tidal captures can hap- 
pen in numbers there, the majority of these would actually be col- 
lisions. This is consistent with our results for encounters with stars 
of this mass. 




10 100 
radius, arcsec 

Figure 9. The observed radial distribution of the X-ray sources in the 
bulge of M3 1 , compared with the expected contributions of different sub- 
populations of low-mass X-ray binaries: primordial (red, the broadest of 
the three histograms), binaries in globular clusters (green, with zero source 
density in the center) and binaries formed through dynamical interactions 
in the inner bulge of M31 (blue, the most peaked distribution). The total 
numbers of sources are summarized in table|6] 



Type 


Number 


Background sources 


29 


Primordial LMXBs 


64 


LMXBs in globular clusters 


21 


LMXBs dynamically formed in the bulge 


21 



6 CONCLUSIONS 

We have studied the spatial distribution of the luminous X-ray point 
sources (L,- > 10"' erg s"') in the bulge of M31. We show that there 



Table 6. Numbers of X-ray sources of different origin in the bulge of M31, 
r < 12 arcmin, Lx > lO'* erg/s 



is a significant increase in the specific frequency of sources, per unit 
stellar mass, in the inner « 1 arcmin. This behaviour is similar, al- 
though smaller in the magnitude, to that observed in globular clus- 
ters. The radial distribution of the surplus sources follows the pi 
profile. All these suggest that the surplus sources are dynamically 
created in stellar encounters in the high stellar density environment 
of the inner bulge of M31. This is further confirmed by the pecu- 
liarity of their luminosity distribution, which res embles that of the 
globu lar cluster sources in M3 1 and our Galaxy jVoss & Gilfanov! 
|2007|) . 

It has long been known that dynamical interactions are respon- 
sible for the relatively large number of X-ray sources observed in 
globular clusters, but this is the first evidence of the dynamical for- 
mation of LMXBs in the vicinity of a galactic center. The stellar 
velocities in bulges are higher than in globular clusters by a factor 
of ~ 5-10. We therefore performed a detailed study of the velocity 
dependence of the three main dynamical processes leading to the 
formation of LMXBs: tidal captures of a compact objects by main 
sequence stars, collisions between evolved stars and compact ob- 
jects and the exchange of a compact object into an already existing 
binary. Another major factor affecting the overal encounter rates 
and the numbers of dynamically formed LMXBs is the high effi- 
ciency of the mass segregation in globular clusters, which modifies 
significantly the spatial distributions of objects of different mass 
and affects the present day mass function in different parts of a 
globular cluster. In addition, due to the relative shallowness of the 
potential well, the populations of compact objects are significantly 
depleted in globular clusters. 

We found that while exchange reactions are potentially the 
dominant formation channel in globular clusters (although stellar 
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collisions might decrease the importance of this channel signifi- 
cantly), this process is relatively unimportant in M31. Similarly, 
tidal captures of NSs by main sequence stars of mass > O.SMq are 
important in globular clusters, but not in M31. Instead the main 
formation channel is tidal captures of compact objects by low mass 
(< 0.3Mq) stars, with some contribution from collisions between 
red giants and compact objects. While the geometrical collision rate 
is high enough to explain the total number of sources from the latter 
channel, the majority of the collisions are unlikely to lead to the for- 
mation of a binary system, as the binding energy of the envelopes 
of most RGB/AGB stars is too low to capture a compact object in 
a high velocity environment. We conclude that the majority of the 
sources in M31 are short-period binaries, and in contrast to globu- 
lar clusters many of them have BH accretors. We note that the BH 
binaries with a very low mass companion may become persistent 
X-ray sources after the donor star is depleted below, Mj < 0. 15Mo, 
due to small size of the accretion disc and the positive dependence 
of the gravitational breaking rate on the mass of the primary. We 
also predict for M31 a large number of faint transients, similar to 
the accreting msec pulsars in our Galaxy. Overall, we have been 
able to explain the spatial distribution and absolute numbers of sur- 
plus sources detected in the inner bulge of M3 1 as a well as the 
statistics of LMXBs in the metal rich globular clusters. However, 
we overpredict by a factor of ~ 1.5 the population of LMXBs in the 
metal poor clusters. 

Finally, the sub-populations of low-mass X-ray binaries in the 
bulge of M3 1 are summarized in Fig|9]and Table[6] 
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